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Abstract. We compare the role of the three-site terms added to the t-J models in the St/ (2) and 
the Ising cases in the extremely low doping regime, i.e. when a single hole added to the strongly 
interacting half-filled system becomes a polaron. We show that in the realistic Ising case the three- 
site terms play a vital role in the polaron movement and should never be neglected unlike in the 
SU(2) case. 
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INTRODUCTION 

The Hubbard model is regarded as the model which can describe many of the spectacu- 
lar phenomena observed in transition metal oxides (TMO) - such as e.g. Mott insulating 
ground states, antiferromagnetism, and perhaps even high temperature superconductiv- 
ity [1]. Usually, however, it is the t-J model which is used to explain these phenomena 
[1]. This model: (i) is easier to solve than the Hubbard model [2], (ii) can be easily re- 
duced to the Heisenberg model for the half-filled doping and hence naturally predicts the 
antiferromagnetic ground state of undoped TMO, (Hi) is believed to give qualitatively 
similar predictions as the Hubbard model for it merely follows from a perturbative ex- 
pansion of the Hubbard model in the physical regime of strong interactions [3]. 

A priori, one can have some doubts concerning the stated above 1 : 1 correspondence 
between these two models. Actually, in a rigorous perturbative expansion of the Hubbard 
model one obtains the so-called three-site terms, in addition to the kinetic and magnetic 
terms in the t-J model [3]. These terms are often neglected in the small doping regime 
(close to the Mott insulating ground state) when they are much smaller than the two t-J 
model terms and only yield a small longer range hopping. Indeed, it was shown in Ref. 
[4] that the inclusion of the three-site terms in this standard SU (2) symmetric case of 
the t-J models does not change qualitatively the low-energy physics of the system. 

However, a different situation can arise when the Ising limit of the t-J model is taken 
and the SU (2) symmetry is broken, which is quite often done for computational purposes 
in e.g. more elaborate t-J models [5]. Could the three-site terms play an important role 
and entirely change the solutions in the Ising limit? In this paper we answer this question 
by comparing the role of the three-site terms added to the t-J models in the SU (2) and the 
Ising cases, cf. [6]. We perform calculations in the extremely low doping regime namely 
for the case of a single hole added to a half-filled ordered ground state of the models, i.e. 
when a single hole couples to the excitations of the ordered state in a polaronic way. 



t-J MODELS WITH THREE-SITE TERMS 



In order to see how the three-site terms arise in the perturbation theory and why they are 
often neglected we introduce (a generalized version of) the Hubbard model: 

H = ~ l £ z fj a ' c L c j<y' + U^n il n i2 , (1) 
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where we defined n ia = c\ a Ci G , a = 1,2 - a quantum number with two values (not 
necessarily a spin quantum number), and zf? is a hopping element between state a on 
site i and state <j' on site j. Hence, this model allows for e.g. a-dependent values of the 
hopping elements. As already mentioned, usually in realistic systems the hopping t is 
much smaller than the on-site interaction U so that we can safely perform a perturbative 
expansion in the kinetic terms [3]. We obtain the effective Hamiltonian: 
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where cJ CT = c] c {\ - n^) and J = 4t 2 /U. 

The first two terms in Eq. (3) form a generalized version of the t-J model (H t -j) 
whereas the last one is the three-site term (H^). The physical interpretation of these 
terms is the following: (i) the first term describes the hopping of electrons with no double 
occupancies allowed and this term contributes to the total energy of the system as oc t8, 
where 8 is the number of doped holes (by definition we have no holes in the half-filled 
case), (ii) the second term oc 7(1 — 8) describes the electron which makes a virtual 
hopping t from site j to an occupied site i which costs energy U and then comes back to 
site j, (Hi) the third term °= JS describes the situation when the electron does not return 
to the same site j but goes to the site / provided this site is unoccupied (hence in the 
half-filled case this term vanishes). Since in the low doping regime 8 1 and J < t, then 
naturally the last term has the smallest contribution to the total energy of the system and 
this is the reason why (as explained in the introduction) it is often neglected. 

Let us also note the reasons for introducing the generalized version of the Hubbard and 
t-J models. Firstly, it enables us to treat the Ising and the SU (2) cases on equal footing 
(see below). Secondly, it is possible to obtain the Ising limit from this model in such a 
way that it coincides with the model describing the correlated electrons in the systems 
with partial t2 g orbital degeneracy (with active yz, zx orbitals in the ferromagnetic plane) 
as e.g. in the d l system of Sr2V04 [7]. This means that this specific Ising limit would 
not only correspond to a standard Ising spin approximation but would also constitute a 
realistic physical model [6]. 




FIGURE 1. Spectral density A(k,o) of the model Eq. (5) with J = 0.4? along the T - M 
(from bottom to top) direction of the 2D Brillouin zone without (with) the first two terms of the 
three-site terms Eq. (6) shown on the left (right) panel, respectively. 



POLARON IN THE SU(2) CASE 

In the SU(2) symmetric case the quantum numbers o =t,| describe physical spins. 
Then 

<4 > 

where the hopping is restricted to the nearest neighbour sites (ij). Substituting Eq. (4) 
to Eq. (3) we obtain for the t-J part 
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where Sf = cj|.c,-j_, SJ~ = c^c^, Sj = \{n^ — n^). The three-site terms part reads 
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where we sum over all possible configurations of three adjacent sites {j'ij} (with i 
denoting the middle site). 

In the half-filled case it is easily seen that the ground state of the model is a quantum 
antiferromagnet (AF). The calculation of the ground state of the model with one addi- 
tional hole in the AF is more complicated and one needs to use the self-consistent Born 
approximation (SCBA) [8]. The hole spectral density A(k, (6) obtained in this way for 
the two-dimensional (2D) lattice is shown in Fig. 1(a) [Fig. 1(b)] without [with] three- 
site terms included, respectively. In both figures we see a finite though strongly renor- 
malized dispersion of the lowest peak. This sugests a formation of a mobile polaron as a 
hole moves by dressing up with spin excitations and acquires a large but finite effective 
mass. The inclusion of the three-site terms Eq. (6) merely strongly increases the spectral 
weight of the incoherent part for particular k values in the Brillouin zone, cf. Fig. 1(b). 
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FIGURE 2. Spectral density A(k,») of the model Eq. (8) with 7 = 0.4/ along the T - M 
direction (from bottom to top) of the 2D Brillouin zone without (with) the first two terms of the 
three-site terms Eq. (9) shown on the left (right) panel, respectively. 



POLARON IN THE ISING CASE 

To obtain the Ising limit from the model Eq. (3) which at the same time serves as the 
realistic model for correlated electrons with ti g degeneracy (see last paragraph of the 
second section) we assume a = yz, zx and 
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where the hopping is between the nearest neighbour sites along the x and y directions 
depending on the value of o [6]. Substituting Eq. (7) to Eq. (3) we obtain for the t-J part 
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where Tf = \{n,i a — The three-site terms part reads 
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where the three-site hop can be either straight (along f or y direction) or along the corner 
[e.g. first hop along the f direction and then along the y one, cf. the last term of Eq. (9)]. 

In the half-filled case the ground state of the model is a classical state with alternating 
orbitals occupied on each site (Ising AF in the spin language). The ground state of the 
model with one additional hole in the Ising AF can again be calculated using the SCBA. 
The hole spectral densities A(k, (6) for the 2D case are shown in Fig. 2 with and without 
three-site terms included. We see that in the case without the three-site terms the polaron 
is immobile [8] whereas adding the three-site terms yields a small dispersion. 



CONCLUSIONS AND FINAL DISCUSSION 



In conclusion, we studied a problem of a single hole doped into the half-filled ground 
state of two different cases of the t-J model. On one hand, in the SU (2) case we noticed 
a quite big quantitative difference in the polaron behaviour between the models with and 
without the three-site terms. However, the polaron was always mobile and the inclusion 
of the three-site terms did not change qualitatively the low-energy physics of the system, 
cf. [4]. On the other hand, in the Ising case the polaron became mobile only after 
including the three-site terms in the model. 

This means that for a t-J model in the Ising limit to become more realistic one should 
include the three-site terms even in the extremely low doping regime. In particular, one 
should also include these terms for the models describing the t2 g systems such as e.g. 
Sr2V04 where the Ising limit arises not as an approximation to the SU (2) case but as a 
genuine t-J model. 

As a postscriptum one can discuss the reason why the three-site terms play any role 
in the Ising limit of the t-J model despite the fact that they are very small in the low 
doping regime (as stated in the introduction and explictly shown in the second section). 
We suggest the following answer: if there is a strong competition between two processes 
in the system (e.g. the magnetic interaction and the hole hopping in the t-J model) then 
the neglected third term (e.g. the three-site hopping) can decide about the ground state 
of the system. This is because it is the difference between these two larger terms in the 
Hamiltonian which matters and with which the magnitude of the neglected term should 
be compared. 
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